The Dynamics of Efficiency: A Simple Model 
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We propose a simple model that describes the dynamics of efficiencies of competing agents. 
Agents communicate leading to increase of efficiencies of underachievers, and an efficiency of each 
agent can increase or decrease irrespectively of other agents. When the rate of deleterious changes 
exceeds a certain threshold, the economy falls into a stagnant phase. In the opposite situation, 
the economy improves with asymptotically constant rate and the efficiency distribution has a finite 
width. The leading algebraic corrections to the asymptotic growth rate are also computed. 

PACS numbers: 05.40.-a, 05.70.Ln, 87.23. Ge 



43 

cd 

CD 



CO 



I 

C 

O 
O 



> 

(N 
(N 

o 
o 

cd 



-a 
c 
o 

CD 



X 



Non-equilibrium statistical mechanics is being increas- 
ingly applied to diverse fields outside physics, ranging 
from biology and computer science to finance and social 
science |l]-|!J. Indeed, the framework of non-equilibrium 
statistical mechanics is ideally suited for describing sys- 
tems composed of many units that interact according to 
simple rules and exhibit a complex large-scale behavior. 
Thus, the important task is to construct simple stochastic 
models incorporating basic characteristics of the dynam- 
ics of systems under study which can then be analyzed 
by employing existing tools of non-equilibrium statistical 
mechanics. The hope is that these models can reproduce 
essential features of the original systems and can help to 
formulate relevant questions and enhance understanding 
of the dynamics of these systems. 

In this paper, we propose a simple model that mimics 
the dynamics of efficiencies of competing agents. These 
agents could be airlines, travel agencies, insurance com- 
panies, etc. In today's competing global economy, the 
performance of a company is continuously judged in the 
market and the index of performance depends on how 
efficient the company is. Rather than trying to incorpo- 
rate all details of performances of competing agents, we 
choose a model that accounts for the dynamics of effi- 
ciency in the simplest form. We represent the efficiency 
of each agent by a single nonnegative number. The effi- 
ciency of every agent can, independent of other agents, 
increase or decrease stochastically by a certain amount 
which we set equal to unity. In addition, the agents in- 
teract with each other which is the fundamental driving 
mechanism for economy. We assume that the interaction 
equates the efficiencies of underachievers to the efficien- 
cies of better performing agents. 

The efficiency model formalizing the above features is 
defined as follows. Let hi(i) is the efficiency of agent i at 
time t. Efficiencies hi's are non- negative integer numbers 
which evolve stochastically. Specifically, in an infinitesi- 
mal time interval At, every hi(t) can change as follows: 



efficiency to that of a better performing agent in 
order to stay competitive. 

(ii) hi(t) — > hi{t) + 1 with probability pAt. This in- 
corporates the fact that each agent can increase 
its efficiency, say due to innovations, irrespective of 
other agents. 

(iii) hi(t) -> hi{t) - 1 with probability q0(hi(t))At, 
where 0(x) is the Heaviside step function. This 
corresponds to the fact that each agent can loose 
its efficiency due to unforeseen problems such as 
labour strikes. Note, however, that since hi{t) > 0, 
this move can occur only when hi(t) > 1. 

(iv) With probability 1 - [1 + p + q9(hi(t))]At, the effi- 
ciency hiit) remains unchanged. 

This efficiency model exhibits rich phenomenology. In 
particular, the system undergoes a delocalization phase 
transition as the parameters p and q are varied. There 
exists a critical line p c (q) in the (p, q) plane such that 
for p > p c (q), the average efficiency increases linearly 
with time, (h) ~ vt for large t. We shall determine ex- 
actly the rate v(j>, q) of the average efficiency growth. For 
P < Pc(q), the economy is stagnant, i.e., the efficiency dis- 
tribution becomes stationary in the large time limit. This 
delocalization (or depinning) phase transition is dynam- 
ical in nature and is different from the depinning transi- 
tions found in equilibrium systems with quenched disor- 
der. Similar delocalization transitions have recently been 
found in a variety of non-equilibrium processes {| ||] . 

Let P(h,t) denotes the fraction of agents with effi- 
ciency h at time t. One can easily write down the evo- 
lution equation for P(h, t) by counting all possible gain 
and loss terms. For h > 1, this equation reads 



dP(h,t) 
dt 



(i) hi(t) 



max[hi(t),hj(t)] with probability At, 



where the agent j is chosen randomly. This move 
is due to the fact that each agent tries to equal its 



P(M) E P(h',t)-(p + q)P(h,t) 

h' = h+l 

qP(h + l,t)+pP(h-l,t) 

h-l 

P(h,t)J2P(h',t). (1) 

h'=0 
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In writing Eq. (|IJ), we have used the fact that when the 
total number of agents diverges, the joint probability dis- 
tribution P(h,h',t) of two agents having efficiencies h 
and b! factorises, P(h,h',t) = P(h,t)P(ti ,t), and the 
mean-field theory becomes exact. 

It proves convenient to consider the cumulative dis- 
tribution, F(h,t) = J2 h '>h p ( h '^)- From E q- (0)> we 

immediately derive the evolution equation for F(h,t), 



dF(h,t) 
dt 



F 2 (h,t) + (l-p-q)F(h,t) 
qF(h + l,t)+pF(h-l,t). 



(2) 



Note that this equation is valid for all h > 1 and by the 
probability sum rule we have F (0, t) = 1 for arbitrary t. 
Also, F(h, t) — > as h — » oo for all t. 

Equation (Q) is a nonlinear difference-differential equa- 
tion and is, in general, hard to solve exactly. Fortunately, 
many asymptotic properties can be derived analytically 
without solving Eq. (||). First we note that F(h,t) ap- 
proaches a traveling wave form as it follows e.g. from 
direct numerical integration of Eq. (^). Thus, we seek a 
solution of the form F(h,t) = f(h — vt). By inserting it 
into Eq. (0) we find that f(x) satisfies 



v^ = -f(x) + (l-p-q)f(x) 
ax 

+ qf(x+l)+pf(x-l), 



(3) 



which should be solved subject to the boundary condi- 
tions /(— oo) = 1 and /(oo) = 0. To determine v, we 
linearize Eq. (|jj) in the tail region, x — > oo. The re- 
sulting linear equation admits an exponential solution, 
f(x) ~ exp(— Ax). By inserting this asymptotics into the 
linearized version of Eq. ([|) we find that the growth rate 
v(X) is related to the decay exponent A via 



v(X) 



l-p 



qe 



pe 



A 



(4) 



Thus we have a family of eigenvalues parameterized by A. 
According to a general selection principle which applies 
to a wide class of nonlinear equations ||[l0]], only one 
specific rate out of this family of possible v's is selected. 
Usually, the minimum rate is selected. For sufficiently 
steep initial conditions, the minimum rate is universal, 
while extended initial conditions might affect the magni- 
tude of the admissible minimum rate. 

The function v(X) in Eq. (Q) has a unique minimum at 
A = A* given by the solution of = 0, or 



1 - p - q + qe A + pe 



,A* 



A* 



-qe 



,-A* 



-A* 



pe 



(5) 



The corresponding minimum rate v m i n (p,q) = v(X*) is 
given by Eq. (0), or 



-qe 



-A* 



pe 



(6) 



as it follows from (|5|). An analysis of Eqs. (Q)-(||) shows 
that there exists a critical line p c (q) in the (p, q) plane, 



p c (q) 



'l + q-2y/q forq>l, 
for q < 1, 



(7) 



such that v(X) > for all A > as long as p > p c (q)- 
For a fixed q, as p — > p c (q) from above, i> m i n — > and for 
< p < p c (q), the curve v(X) crosses zero at A = Ai and 
A = A2 with A2 > Ai. When Ai < A < A 2 , v(X) becomes 
negative. This tells that there might be no traveling wave 
solution for < p < p c {q) and we anticipate that the effi- 
ciency distribution should become stationary. Note that 
for q < 1, p c (q) = and this regime does not occur. 




FIG. 1. The thick line represents the critical locus, 
Pc(q) = 1 + q — 2y/q, for the mean-field theory. The +'s 
indicate numerically obtained critical points in one dimen- 
sion. For sharply decaying initial conditions, the economy is 
developing when p > p c (q) and stagnant when p < p c (q). 

With the above picture in mind, we now discuss the 
selection principle more carefully. Consider an exponen- 
tially decaying initial condition, F(h, 0) ~ e~ ah with 
a > 0. When p > p c (q), the rate is positive for all 
A > and v(X) has a unique minimum at A — A*. Ap- 
plying the selection principle we find that for sufficiently 
steep initial conditions, a > X* , the selected growth rate 
is Wmin = w (^*)- Consider now sufficiently extended ini- 
tial conditions, a < A*. In this case, f(x) must decay at 
most as e~ ax and therefore the growth rate is selected 
among v(X), Eq. (||), with A < a. The selection principle 
now implies that the selected rate is v = v(a). 

When p < p c {q), we must separately consider two 
cases: q > 1 and q < 1. For q > 1, v(X) as given by 
Eq. ([|) becomes negative in the region Ai < A < A 2 . We 
find that for all a < X\, the system still admits a travel- 
ing wave solution and the selected rate is v = v{a). How- 
ever, for a > Ai, the system no longer admits a traveling 
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wave solution. Instead, the distribution F(h, t) reaches a 
stationary limit Poo{h) as t — > oo. By putting the time 
derivative equal to zero on the left-hand side of Eq. (0), 
we find that the stationary efficiency distribution decays 
exponentially, F^h) ~ e _Al/l , with 



MO, q) = ln 



-1 + p + q + a/(1 — p — q) 2 — 4pq 



(8) 



Note that [i(p, q) is real below the critical line, i.e., when 
q > 1 and p < p c {q)- Interestingly, q) remains finite 
on the critical line p = p c (q)- From Eqs. (||) and (Q) we 
find that fj, c (q) = n(p c (q),q) reads 



/i c (q) = ln 



(9) 



For q < I, Pc{q) = 0. When p — > 0, we have u m i n — * 
and A* — * 00. The divergence of the decay exponent A* 
indicates that when p = and q < I, the system still 
admits a traveling wave solution and the selected rate is 
v = v(a) if we start with an exponentially decaying ini- 
tial condition, F(h,0) ~ &- ah . Of course, for compact 
initial conditions (i.e., when F(h,Q) — for sufficiently 
large h), the efficiency distribution becomes stationary 
in the long time limit. We have verified all the above 
assertions via direct numerical integration of Eq. (^) . 

Although one cannot provide explicit expressions for 
the growth rate in the developing phase, near the critical 
line p c (q) the growth rate considerably simplifies. First 
we note that on general scaling grounds one would guess 
that above the critical line, the growth rate t> m in(p, q) 
should be a function of p — p c with critical behavior 



(p ~ Pc) 
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(10) 



The actual behavior is found by a straightforward anal- 
ysis of Eqs. (@)-(g) to yield 

p-Pc(q) 



ln(l/p) 
1-g 



for q > 1 , 
for q = 1, 
for < q < 1 , 



(11) 



ln(l/p) 

where /i c (<z) is given by Eq. (ph. Equation ( [Tl"| ) implies 
that the mobility exponent (3 in the scaling relation jlf| ) 
is equal to 1,1/2 and for q > 1 , q — 1 and q < I , re- 
spectively. In the last situation (q < 1 and p — * 0), the 
growth rate still approaches to zero but it occurs in an 
extremely slow inverse logarithmic fashion. 

The relaxation of the growth rate v(t) towards its 
asymptotic value v m i n exhibits an interesting algebraic 
behavior. Specifically, the leading correction is propor- 
tional to the next is of order £~ 3 / 2 , etc. Similar 
i _1 correction was first derived by Bramson in the con- 
text of a reaction-diffusion equation [ fio| , and was sub- 
sequently re-derived and generalized by a number of au- 
thors The next correction was recently derived 



by Ubert and van Saarloos p3). In contrast to Refs. 
]To| p^[ , we consider the difference-differential equation. 
Fortunately, the techniques JlO) |l3|] still apply (compare 
HQ), so we do not detail the derivation. Following for 
instance an approach of Ref . [|l3| , one finds 

v(t) = v min - JL t -i + At- 3 / 2 + 0{r 2 ), (12) 

with A = 37r 1 / 2 {2{qe- x * +pe^)}- 1 / 2 (\*)- 2 . The explic- 
itly displayed terms are universal - they do not depend on 
initial condition as long as it is steep enough [i.e., it falls 
off faster than e _A x ]. The following terms in Eq. jf^ ) 
starting from 0(i~ 2 ) correction are non-universal. Thus 
not only any sufficiently steep initial profile relaxes to a 
unique profile, the approach to that profile occurs along 
(asymptotically) unique trajectory. Note also that the 
very slow i _1 relaxation of the growth rate leads to a 
logarithmic correction to the average efficiency, 



(h) 



2A" 



lnt + 0(l). 



(13) 



Thus, we have a derealization transition across the 
critical line p c (q) in the (p, q) plane for sharply decaying 
initial conditions. For such initial conditions, as long as 
p > p c (q), the economy is in the developing phase with 
the average efficiency increasing as (h) ~ v m i n t, where 
the rate u m i n given by Eqs. (||) and (||). For p < p c (q) 
with q > 1, the system is localized and (h) approaches 
a time-independent constant in the long time limit. For 
p = and q < I , the economy is in the developing phase 
for unbounded initial efficiency distributions, with the 
rate of growth explicitly dependent on the initial con- 
dition. For economically more relevant compact initial 
conditions, the regime p — and q < I belongs to the 
stagnant phase. These results are summarized in the 
phase diagram in Fig. 1. 

The mean-field version of the efficiency model is nat- 
ural in the interconnected modern economy. In econ- 
omy with limited communication, however, the efficiency 
model in a low dimensional space rather than in the fully 
connected graph might be more appropriate. In this case, 
agents are placed on a finite dimensional lattice. The mi- 
croscopic dynamical steps (i)-(iv) remain the same ex- 
cept that in move (i), the agent j is chosen to be one 
of the nearest neighbors of i. Unlike the mean- field the- 
ory, the correlations between hi's at different sites remain 
nonzero in finite dimensions even in the thermodynamic 
limit. We have studied this model numerically in one di- 
mension. The results are shown for lattice size L = 1 000 
(we verified that for such large systems, the finite size 
effect is insignificant). Once again, there is a derealiza- 
tion transition in the (p, q) plane across a critical line as 
shown in Fig. 1. The efficiency distribution P{h,t) at 
different times in both phases is presented in Fig. 2. 
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FIG. 2. The distribution P(h, i) at times t = 200, 600 and 
1000 in the moving phase for p — 3, q — 4. The inset shows 
the distribution at the same times in the localized phase for 
p = 1, q = 4. For q = 4, the critical point is p c (q) ~ 1.7. 



models in (l + l)-dimensions belong to the directed perco- 
lation (DP) universality class (see e.g. Ref. [T^j). Other 
similar growth models exhibit phase transitions that do 
not belong to the DP universality class It remains 

an open question whether the phase transition in the ef- 
ficiency model in (1 + l)-dimensions belong to the DP 
universality class. 

In summary, we have investigated a simple model of the 
dynamics of efficiencies of competing agents. The model 
takes into account stochastic increase and decrease of the 
efficiency of every agent, independent of other agents, 
and interaction between the agents which equates the ef- 
ficiencies of underachievers to that of better performing 
agents. We have shown that the model displays a phase 
transition from stagnant to growing economy. 

One of us (PLK) acknowledges support from NSF (grant 
DMR9978902) and ARO (grant DAAD19-99-1-0173). 



We now stress important differences between mean- 
field and finite-dimensional situations. In the former case 
the nature of the two phases depends on the steepness pa- 
rameter a, while in one dimension the nature of the final 
state is independent of a. For example, in the develop- 
ing phase the system always has a traveling wave solution 
with a rate that depends on p and q but does not depend 
on a. We have tested this fact numerically for several 
values of a. This result is rather counter intuitive as 
it suggests that correlations seem to restore universality 
that mean-field theory lacks. Another important distinc- 
tion is a very different behavior of the width of the ef- 
ficiency distribution in the developing phase. Indeed, in 
mean- field the width is constant, while in one dimension 
it increases with time [see Fig. 2]. Moreover, the width 
increases as a power law, to = ^/{h 2 ) - (h) 2 ~ for 
large t, with (3 w 0.31 in (1 + 1) dimensions. 

In (d+ l)-dimensions, one can interpret the efficiency 
hi(t) as the height of a surface growing on a d-dimensional 
substrate. In this language, our efficiency model rep- 
resents a continuous time polynuclear growth (PNG) 
model with additional adsorption and desorption rates 
flBf . The continuous PNG model without desorption has 
been studied within mean-field theory ||l6[| and was found 
to be always in the moving phase as expected. From 
the general analogy to PNG models, we expect that the 
moving phase in the efficiency model corresponds to a 
growing interface belonging to the Kardar-Parisi-Zhang 
(KPZ) universality class fis|| . The numerically obtained 
width exponent (3 sa 0.31 in (1 + l)-dimensions is con- 
sistent with the KPZ prediction (3 — 1/3. It would be 
interesting to determine the universality class of the de- 
localization transition. Phase transitions in several PNG 
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